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Abstract: We revisit bosonization of non-relativistic fermions in one space 
dimension. Our motivation is the recent work on bubbling half-BPS geometries 
by Lin, Lunin and Maldacena (hep-th/0409174). After reviewing earlier work on 
exact bosonization in terms of a noncommutative theory, we derive an action for 
the collective field which lives on the droplet boundaries in the classical limit. 
Our action is manifestly invariant under time-dependent reparametrizations of the 
boundary. We show that, in an appropriate gauge, the classical collective field 
equations imply that each point on the boundary satisfies Hamilton's equations 
for a classical particle in the appropriate potential. For the harmonic oscillator 
potential, a straightforward quantization of this action can be carried out exactly for 
any boundary profile. For a finite number of fermions, the quantum collective field 
theory does not reproduce the results of the exact noncommutative bosonization, 
while the latter are in complete agreement with the results computed directly in 
the fermi theory. 
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1. Introduction 



The connection between free non-relativistic fermions and string theory in 2- 
dimensions is known since early nineties 1 [ffL |2|, |3|, [|, |^, [], [7], ^ [| . Recent studies 

have shown that free non- 



14, 15, 16, 17, 18, 19 



1L 22, E3L 24, 25 



relativistic fermions also appear in other situations in string theory, typically in 
sectors which have high enough supersymmetry. The requirement of sufficient 
amount of supersymmetry is understandable since 2-dimensional free fermions form 
an integrable system. The correspondence with fermions is usually more easy to 



see in the dual field theory. Remarkably, in [T3|] the authors found a class of 
1/2— BPS solutions of supergravity equations in one-to-one correspondence with 

1 For an older review of the subject see ftc[ |; more recent review and references can be found 
in 0,R3. 
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classical configurations of free non-relativistic fermions in 2-dimensions. While 
providing yet another example of AdS/CFT correspondence 2 , it opens up the 
interesting possibility of learning something about the nature of quantum gravity 



and string theory from free fermions |28], |23[. Since in [13|] fermions make contact 
with geometry via the bosonized theory which describes their collective motion, it 
is essential to understand all aspects of the bosonized theory in order to be able 
to use the full power of free fermions. This provides the main motivation for the 
present work. 

The organization of this paper is as follows. In Section 2 we review (see also 



]) the works |29], [3C], |31| in which an exact bosonization of free non-relativistic 
fermions in 2-dimensions has been developed in terms of a noncommutative theory 
using Wigner phase space density. In Section 3 we discuss the classical limit of 
this bosonized theory. In the classical limit, a generic configuration consists of 
droplets of fermi fluid on the phase plane. The dynamics, which is associated to 
the collective motion of the droplets, manifests itself in their changing boundaries. 
The action for this collective motion has a built in symmetry under arbitrary time- 
dependent reparametrizations of the droplet boundaries. This gauge symmetry 
reflects the fact that the motion of the fluid along the boundaries of the droplets 
is unphysical because of the indistinguishability of fermions. In this section, we 
derive a manifestly gauge-invariant classical action which describes the boundary 
dynamics of the droplets for arbitrary boundary profiles. Issues related to gauge 
fixing for different droplet boundary profiles are also discussed in this section. In 
an appropriate gauge, we show that classically each point on the boundary simply 
follows Hamilton's equations of motion. Quantization of the collective field theory 
is carried out in Section 4. For the harmonic oscillator potential, quantization can 
be carried out exactly for any droplet boundary profile. We find that for a finite 
number of fermions, the spectrum of the collective quantized theory does not agree 
with the exact spectrum at large energies. Furthermore, the phase space density 
fails to reproduce the precise details of the exact result. We summarize our results 
and end with some concluding remarks in Section 5. 

2. Review of exact bosonization 

Consider N free non-relativistic fermions moving in one space dimension in a 
potential. The fermion wavefunctions satisfy the Schroedinger equation 

ih d t ip{x,t) = Hip{x,t), (2.1) 

where the single-particle Hamiltonian H is given by 

H= 1 -(-h 2 dl + V(x)). (2.2) 



2 For a review see [E7j. 
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Moreover, the fixed fermion number constraint is, 

r+oo 

dx 

m=l 



r+oo N 

N= dx Y,^m( X ^m{x,t). (2.3) 
J —oo . i 



Here ip m (x,t) (m = 1,2,3- ■ •) form a complete orthonormal set of single-particle 
wavefunctions, i.e. S mn = J_££ dx i[)^(x,t)i/; n (x,t). 

2.1 Bosonization in terms of Wigner density 

The bosonization carried out in [^, [30], |3lJ uses the Wigner phase space density 



as its basic building block. In terms of the fermion wavefunctions, it is defined by 
the expression 

r+oo N 

u(p,q,t)= dxe-^ n *E^(q-x/2,t)iJ> m (q+x/2,t). (2.4) 

J -°° m=l 

A fundamental property of the Wigner density, one that captures the 
underlying fermionic nature of the degrees of freedom in the bosonized version, 
is that it satisfies a nonlinear constraint. This constraint, which can be elegantly 
written using the non-commutative star product 3 in the phase plane, is 

u*u(p,q,t) = u(p,q,t), (2.5) 

where the star product is defined in the usual way, 



Mi * u 2 (p,q,t) 



e f(d qi a P2 d i2 d Pi) Ul ( Puqu t)u 2 (p2,q2,t) 



■ (2-6) 

91=92=9, Pi=P2=P 



A quick way of deriving the constraint is to first construct the bilocal fermion 
bilinear J2m=i ^mfoi ttymiy, t) = <p(x, y, t). As a consequence of the orthonormality 
of the fermion wavefunctions, this bilocal function satisfies the constraint 
f_™ dz (f)(x, z,t)(f)(z,y,t) = (f)(x,y,t). The definition of u, eqn( [2.4| ), which can 
be easily rewritten in terms of the bilocal function (p(x,y,t), and some simple 
algebraic manipulations then lead to the constraint ( |2.5| ) on u. 
In addition to the above constraint, u satisfies the condition 

P-7) 

which is just a restatement, in terms of the Wigner density u, of the fact that the 
total number of fermions is N. Finally, an equation of motion can be derived for u 
using the Schroedinger equation satisfied by the fermion wavefunctions. One gets, 

dtu(p, q, t) = {h, u}*(p, q, t). (2.8) 



3 The compact star product notation was not used in references |^9|, The expressions 

given there are however the same with star products written out in long-hand. 
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Here h — |(p 2 + V(q)) is the classical single-particle Hamiltonian. The bracket on 
the right-hand side is Moyal's generalization of the Poisson bracket involving the 
star product, namely = -^{h * u — u* h). In the limit % — > the Moyal 

bracket goes over to the Poisson bracket. 

2.2 The role of symmetry 



Notice that the constraints fl2.5|) and (|2.7| ) are left unchanged by the following 
infinitesimal variation of u, 

5u = {e,u}«, (2.9) 

where e = e(p,q,t). This is, in fact, the most general variation that leaves these 
two constraints unchanged. If the t-dependence of e is such that d t e = {h, e}* , 
then u + 5u also satisfies the equation of motion ( |2.8| ). It is easy to see that two 
such successive transformations satisfy the group composition law. In fact, the 
relevant group is W Ml the quantum generalization of the group of area-preserving 
diffeomorphisms in 2-dimensions, and these transformations move us on a co- 
adjoint orbit of Woo in the configuration space of w's. In particular, changes in 
u corresponding to t-independent e's satisfying {h, e}* = are symmetries. The 
corresponding conserved charges 4 are 

Qn= /tj h n *u(p,q,t), n=0,l,2,--., (2.10) 

where h n = h*h*---*h has n factors. 

The charges Qq and Q\ are familiar. They measure the total number of 
fermions and total energy, respectively. Other charges are less familiar in the 
bosonic version. In terms of fermions, however, they can be easily seen to be sums 
of higher (than linear) powers of individual fermion energies, which are obviously 
conserved in the non-interacting theory that we are considering here 5 . 

2.3 The action of the bosonized theory 



Equations ( |2.4| )-( |2~lf ) constitute an exact bosonization of the fermion problem. In 
[p9| pO] , |3~T| a derivation of this bosonization has been given and an action obtained 
for the Wigner density, whose variation gives rise to the equation of motion 



4 Similar charges first appeared in Q in a system of free fermions in inverted harmonic oscillator 
potential, which is relevant for the c = 1 matrix model. 

5 In the fcrmionic theory, all the charges in any given state can be explicitly written in terms 
of the individual energies of N fermions. One might wonder how the bosonic version of these 



charges in ( [2.10 ) depends on only N independent parameters. The point is that the u's that 



should be used to evaluate these charges must satisfy the constraints (|2.5| ) and (2/7). Once this 
is ensured, it can be seen that the charges Q n can be expressed in terms of only N independent 
parameters. 
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The derivation uses coherent states and the coadjoint orbits of W mi much like the 
coherent states and coadjoint orbits of SU(2) are used in the case of a spin in a 
magnetic field. As in the latter case, one needs to construct a "cap" action, the 
cap being parametrized by time t and an additional variable s (0 < s < 1), such 
that u(p, q, t) = u(p, q,t;s = 1) while u(p, q, t; s — 0) = u(p, q) is a t-independent 
function. One gets the following action: 

S = JdtJ ^ (J ds h 2 u * {d t u, d s u}* -u*h)j (2.11) 

It can be easily verified that the equation of motion (|2.8|) follows from this action 
upon using the variation ( |2.9| ) and the constraints ( |2.5|) and ( |2.7| ). 

A simpler form for the action can be obtained, one that is useful for going over 
to the classical limit, if one rewrites the action in terms of a "reference" density 
u (for example, it could be the density in the fermi vacuum) and the W m group 
element v 6 that is needed to "rotate" uq to u = v * uq * v'. One gets 

where a t = ih dtv * v'. Since at also satisfies the equation 

d t u = {a t ,u}*, (2.13) 

which does not depend on the reference density Uq, we may alternatively use this 
equation to define a t . We will later use the equations ( j2.12j ) and ( |2.13j ) to take the 
classical limit of this bosonized theory. 



2.4 Consistency of bosonization 

We will end this section by arguing that the equations ( |2.5| )-(|2~8|) provide an 
exact bosonization in the sense that solutions to these equations are in one-to-one 
correspondence with the states of the fermion system. Consider first the constraint 



equation ( [2.5|) . General solutions to this equation have been obtained in p3| . For 
the harmonic oscillator potential for which h = \{p 2 + q 2 ) = p, which is the case 
of interest in |13|], the analysis is particularly simple. This is because in this case 



t-independent solutions of (|2.8|), which correspond to energy eigenstates, are radial 
functions in the phase plane. As discussed in |33), there are an infinite number of 
real independent radial solutions to the equation ( |2.5|) . The general form of the 
solution is 



m=0 



G v is unitary, that is v * v' = * v = I. 
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where <fi m (p) = 2(— l) m e 2h L m (Ap), L m being the mth Laguerre polynomial. The 
coefficients take the values 0, 1. The fermion number constraint ( |2.7|) fixes 

the number of nonzero c m 's to be precisely N. These, in fact, identify the filled 
levels, so it is clear that solutions to the equations (f2.5|)-(p~8|) are in one-to-one 
correspondence with the energy eigenstates of the fermion theory. In fact, given a 
time-independent u(p,q) which solves the equations (|2T5|)-(|2"78|), one can uniquely 
reconstruct the filled levels. Since u is essentially a projection operator, finding the 
corresponding fermi state amounts to finding the subspace on which it projects: 

^ e«-W* u{p, X -±«) 4,{q) = ^(x). (2.14) 



For N fermions, this subspace is N- dimensional, so equation (|2.14 ) should have N 



linearly independent solutions. In this way one can obtain the wavefunctions for 
the filled states from the given u. 

We give below the solution for the Wigner density in the fermi vacuum: 



N-l 

u F =Y. 2(-l) m e- 2 ^L m (4p). (2.15) 

m=0 



The expression on the right hand side can be rewritten 7 as an integral over a single 
Laguerre polynomial: 

u F = dye-y^L^M- (2-16) 

J&p/K 



Note that for large N p4| , 

e-» /2 4r-i(l/) = 7r 1/2 y- 3/4 W 1/4 cos(2yJ(N-l)y-3ir/4) + 0(A^ 1/4 ). (2.17) 

We will make use of these expressions when we compare u p given above with the 
one obtained by quantizing the classical limit of the action in (|2.13| ). 

3. Classical limit and collective action 

The easiest way to take the classical limit of the bosonized theory obtained as 
outlined above is through the equations (|2.12|) and fl2.13D , supplemented by the 



constraints ( |2.5| ) and ([2.7]) . To be precise, in the classical limit we will set h — » 
and — > oo while keeping ^.A^ = po fixed. In this limit equation ( |2.5| ) becomes 
u 2 = u, whose standard solutions are filled and empty regions of phase plane. 
Equation ( |2.7| ) fixes the total area of the filled regions, J dpdq u = 2irp . Note also 
that the classical energy of any configuration diverges in this limit as 0(1 /h). 

7 Actually the answer depends on whether N is even or odd. The density falls off to zero at 
infinity only for N odd. It is this answer that we have given in the equation below. 



6 



3.1 Classical motion is area-preserving diffeomorphisms 



Let us now first consider equation (|2.13| ). In the classical limit, the Moyal bracket 
on the right-hand side reduces to the Poisson bracket. Integrating this equation 
over an infinitesimal amount of time St gives 

u(p — St d q a t , q + St d p a t , t + St) « u(p, q, t). 

This equation shows that the classical motion of the fermi fluid is determined 
by area preserving diffeomorphisms in phase space. Here the diffeomorphism is 
generated by the function a t . Physical motion of the fermi fluid is manifested only 
in changing boundaries 8 of the filled regions. Since the point (p, q) in the phase 
plane moves to the point (p — St d q a t , q + St d p at) in an infinitesimal time St, for 
consistency we must have 



d t p = -d q a t , d t q = d p a t . 



(3.1) 



We will use these equations below to get an explicit form for the classical action 
involving only the boundary variables. 

3.2 The boundary action 



Consider the action in ( |2.12| ) in the classical limit. For simplicity we will assume 
that the density u is nonzero only in a single connected region E, centered at 
the origin of the phase plane, as shown in Fig. 1(a). Generalization to several 
disconnected filled regions, as for example in Fig. 1(b), is straightforward. Let 





(a) 



(b) 



Figure 1: (a) A simple configuration with a single connected filled region S in phase 
space. The droplet is centered at the origin in the phase plane. (b) A general 
configuration of several disconnected fluid droplets on the phase plane. 

us parametrize the filled region S by (r, cr) which take values in the unit disc 

8 This is because in the interior of the filled regions, the motion is unphysical since the fermions 
filling the interior are identical. Even on the boundary of a filled region, motion of the fluid along 
the boundary is unphysical for the same reason. 
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< t < 1, < a < 2n. Then the classical action can be written as an action on 
this disc 

S = — !-r f dt f dr f da (d a pd T q - d a qd T p)(a t - h). (3.2) 
2nn J Jo Jo 

For later purposes, it will be more convenient for us to work in polar coordinates 
9 , q = rcosd, p = rsin9. In these coordinates the above action becomes 

S = — ^— [ dt [ dr [ da (d a 9d T p — d a pd T 9)(a t — h), (3.3) 
27m J Jo Jo 



where we have used r 2 /2 = p. We also write down the consistency conditions (|3.1| ) 
in these coordinates: 

d t 6 = -dpOt, d t p = d e a t . (3.4) 

Since physical motion of the fermi fluid is manifested only in changing 
boundary of the filled region E, it should be possible to reexpress this classical 
action in terms of appropriate degrees of freedom which live only on the boundary 
of the disc. To do so, let us introduce the collective field 10 0, which is defined by 
the relation 

pdj = d a <p. (3.5) 

In terms of this variable, 

(d a 9d T p - d a pd T 9) = d a (d T (f) - pd T 9). 

Also, using the consistency conditions (|3.4|) , we get 

d a a t = d t pd a 9 - d t 9d a p 
= d a (d t <j ) -pd t 9). 

With the help of these two relations, and after some algebraic manipulations 
involving partial integrations, the integrands of both the kinetic term involving 
a t as well as that of the hamiltonian piece in action ( |3.3|) can be expressed as total 
derivatives in r. The resulting boundary action (at r = 1) is 



S = — - I dt I da 

2nn 



2 TV 



^d a (j>{d t (j> - P d t 9) - djh 



(3.6) 



where d p h = h 11 . This action, together with the fixed area constraint, / da d a 9p = 
2tcpq, describes the dynamics of the boundary of the filled region S. 

9 For applications to 2-dimensional string theory and the c — 1 matrix model, it would be 
more appropriate to work in hyperbolic coordinates. 

10 The collective field approach to approximate bosonization of non-relativistic fermions was 
first used in |p^| . 

11 For example, for the harmonic oscillator potential, h = + q 2 ) = p and so ft. = p 2 /2, upto 
an irrelevant constant. 
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3.3 Boundary reparametrizations 

Remarkably, the equations of motion arising from the above action for independent 
variations with respect to <fi and 9 turn out to be identical. This equation is 



d a 9d t p - d t 9d a p - d a h = 0. 



(3.7) 



It is not hard to understand the reason for this. The action ( |3.6| ) is invariant 
under t-dependent a-reparametizations, a — > a'{t,a). This gauge invariance of 
S arises because the motion of the fluid along the boundary is unphysical, which 
is due to the indistinguishability of the underlying fermionic degrees of freedom. 
As a consequence of this gauge symmetry, p(t, a) and 9{t, a) provide a redundant 
description of the dynamics of the boundary. The physical description of the 
dynamics requires only "half" the number of variables. The single equation of 
motion (|3.7| ) describes the dynamics of this physical, gauge-invariant degree of 
freedom. 

For a generic potential (|3.7|) is a complicated non-linear partial differential 
equation, not easy to solve. However, for the harmonic oscillator potential it 
simplifies dramatically. In this case h = p, so the equation becomes 



A general solution to this equation can be written down immediately. At points 
on the boundary where neither d a p nor d a 9 vanishes, the solution is 



Here / is any arbitrary periodic function (in 9 — > 9 + 2ir), subject only to the fixed 
constraint / da d a 9p = 2np . 

The above solution is not valid at the points crj where either d a p or d a 9 vanishes. 
At these points (|3~^ ) requires both d a p and d a 9 to vanish together, which is not 
possible 12 . The equation of motion is, however, satisfied everywhere. Therefore, 
at the points where d a p vanishes, we must have d t p = 0. Similarly, at the points 
where d a 9 vanishes, we must have dtO = — 1. To proceed further and completely 
solve the classical dynamics of the equation ( |3.8|) , we need to suitably gauge- fix its 
symmetry under t-dependent cr-reparametizations. This is what we will do now. 

3.4 Gauge-fixing boundary reparametrizations 

We will continue our discussion with the specific example of the harmonic oscillator 
potential. Generalization to other potentials is easy and will be mentioned at the 

12 There is a possibility of both d a p and d a vanishing for self-intersecting curves at the point 
of self-intersection. However, such boundary profiles do not seem meaningful for regions filled 
with fermi fluid. Moreover, a self-intersection of the boundary is a pinching of the filled region, 
which cannot be meaningfully described in the classical limit. 






(3.9) 
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end. We start by choosing a gauge in which p is a fixed t-independent function 
of a. This choice is possible, except at points where d a p vanishes, since p is 
gauge-invariant at these points. However, by equation of motion (|3.8|), p is still t- 
independent at these points. Thus, in this gauge it is consistent with the equation 
of motion to choose p(t, a) = p(a) everywhere. From ( |3.8| ) it follows that we must 
have 9(t, a) = —t + 9(a), except possibly at the points where d a p vanishes. But by 
a gauge choice we can adjust 9 to this solution at these points as well 13 . Hence, 
in this gauge, the complete gauge-fixed solution of (|3.8|) is 14 



p(t,a)=p(a), 9(t,a) = -t + 9(a). (3.10) 

We see that each point on the boundary of the filled region simply rotates in a 
circle around the origin in the phase plane. This is precisely what happens in 
the fermi picture in the classical limit - in the harmonic oscillator potential the 
particles simply rotate in circles whose radii are determined by their energies. 

It is easy to generalize the above argument to other potentials. Let us rewrite 
the equation of motion ( |3.7| ) for a generic potential as follows: 

dj(d t p - d h) - d a p(d t 9 + d p h) = 0. (3.11) 

As above, we gauge-fix p such that its t-dependence is determined by the equation 
dtp—dgh = 0. By ( |3.11| ) this equation continues to be satisfied even at points where 
d a p vanishes and p is gauge- invariant. Then, from (|3.11|) we get d t 9 + d p h = and 
by gauge-fixing 9 we can enforce this equation even at points where d a p vanishes. 
The result is that in this gauge, physical dynamics of the boundary is obtained by 
solving the equations 

d t p = d e h, d t 9 = -d p h. (3.12) 

These are precisely Hamilton's equations for a particle with the hamiltonian h. So 
the points on the boundary of the filled region follow the trajectories described by 
the solutions to these equations. This is exactly as expected from fermions moving 
in the given potential in the classical limit. 

We end this discussion by explaining the physical meaning of the above gauge- 
fixing procedure. At generic points on the boundary of a filled region, the motion of 
the fermi fluid can be arranged to be purely angular or purely radial by adding an 
arbitrary motion of the fluid along the boundary. The symmetry under t-dependent 
cr-reparametizations allows us to do this. At the points cr, where d a p (d a 9) vanishes, 
however, the tangent to the fluid boundary is purely angular (radial), physical 
motion is purely radial (angular) and only angular (radial) motion can be changed 
by gauge changes. We use this latter freedom to adjust, consistent with equation 
of motion, 9(t, a) (p(t, a)) to be a smooth function of t across the points (Xj. 

13 Since p is gauge-invariant at the points where d a p vanishes, we use the freedom of gauge- 
transformations to adjust 9 at these points. 

14 This solution is consistent with (fOl) for values of a for which the latter holds. 
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3.5 Conserved charges and area-preserving diffeomorphisms 

The conserved charges that we found in ( |2.10| ) in the exact bosonized theory have 
a classical limit. In fact, the classical expressions for the charges can be written 
entirely in terms of the boundary variables which appear in the action ( |3.(j| ). We 

get 

Qn = 7TT d ° ~ h m n = 0,1,2,- ••, (3.13) 
2irn Jo 

where d p h n = h n . One can check directly using the classical equation of motion 
( |3.7| ) that these charges are conserved. Note that, unlike the charges in ( |2.10j ), here 



the number of independently conserved charges is infinite. This is consistent with 
the fact that in the classical limit the number of fermions goes to infinity. 

In the case of harmonic oscillator potential, these charges take an especially 
simple form: 

C £ = T% dad c 6f—- 71 = 0,1,2,.-., (3.14) 
Zirh Jo (n + 1) 

In addition, in this case one can give an explicit expression for the other charges 
of the group of area-preserving diffeomorphisms: 

C^r^^^ 1 ^" 1 '- n,m = 0,l,2,..-. (3.15) 
AixTi Jo 

Using the equation of motion ( |3.8|) it can be checked that these satisfy the equation 

d t W* m = i{n-m)W* m . (3.16) 



The "diagonal" charges, W^ n , are conserved. They are, in fact, essentially the 
charges Qf 



4. Quantization of the collective theory 

In this section we will discuss the problem of quantization of the collective theory. 
Throughout this section we will limit our discussion to the harmonic oscillator 
potential. This is because in this case quantization can be done exactly. One can, 
of course, develop a perturbation expansion in small h for more general potentials, 
but we will not pursue that here. 

Our staring point is the classical action (|3.6| ). Canonical quantization requires 
that we first fix a gauge for the symmetry under boundary reparametrizations 
that this action possesses. The gauge that we will fix is slightly different from the 
one we discussed in the previous section. Here we will start by gauge-fixing 9. A 
convenient choice is to set 8(t, a) = —t + 9(cr). By arguments similar to those given 
in the previous section, we again arrive at the gauge- fixed solution ( |3.1C| ). 
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The choice of the precise functional form of 9(a) depends on whether d a 9(a) 
vanishes anywhere or not. Geometrically, the tangent to the boundary at such 
a point is in the radial direction. Physically, at these points 9 is gauge-invariant 
and the gauge is fixed on p, to arrive at the gauge-fixed solution ( |3.10| ). The 
simplest case is that of a boundary which has no such points. More generally, 
however, there may be several points on the boundary at which the tangent is 
radially directed. Clearly, the number of such points has to be even because the 
boundary is a closed curve, unless d 2 9(a) also vanishes at some point, i.e. it is a 
point of inflexion. Additionally, there may be points where still higher derivatives 
of 9(a) also vanish. In the following, we will separately discuss quantization for 
the two cases: (i) tangent not radially directed at any point on the boundary and 
(ii) tangent radially directed at one or more points, some of which may also have 
vanishing higher derivatives of 9(a). 

In the following, for simplicity we will restrict our detailed discussion to a single 
fluid droplet, assuming further that the droplet is centered around the origin in 
phase plane. At the end of this section, we will comment on extension of this 
discussion to more general fluid configurations of the type shown in Fig. 1(b). 

4.1 Boundary profiles with tangent nowhere radially directed 

This is the simplest case. A representative example of this class of boundary profiles 
is shown in Fig. 1(a). By a gauge choice, consistent with equation of motion, in 
this case we may set 

9(t,a) = -t + 9(a), 9(a) = a. (4.1) 

The canonical equal-time commutation relation (actually the Dirac bracket) for (f) 
that follows from the action ( j3.6|) is 

[d a <j>{t, a), d a ,<j)(t, a')] = -2i7rh 2 dJ(a - a'). (4.2) 

Using fl3.5| ) and ( |4.1| ), the commutation relation for p follows: 

[p(t, a),p(t, a')] = -2inh 2 dJ(a - a'). (4.3) 

In the above gauge, p is time-independent. In terms of modes, 

p(t, a) = p(a) = Po + hjr(a m e m ° + al e~ m(7 ) . (4.4) 

m=l 

The constant term is fixed to be po because of the fixed area constraint J da d a 9p = 
2np . The physical degrees of freedom are the complex modes a m (m = 1, 2, • • •). 
Because of (^4.3|) they satisfy the harmonic oscillator commutation relations 

[a m ,o4] = m5 mn . (4.5) 
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The hamiltonian is given by 



H = ^-JdadJ p 2 . (4.6) 
In terms of the modes of p this reads 

H = a ln a m + ttt + "zero — point energy". (4.7) 

m=l 

The second term in the hamiltonian is the ground state energy. In the fermionic 
picture, this is the energy of the fermi ground state. The last term is an infinite 
"zero-point energy". 

The first term in the hamiltonian gives excitation energies. The excited states 
are constructed from an infinite number of decoupled oscillators with frequencies 
1, 2, • • •, just like in free string theory. For low energies, these states are in one-to- 
one correspondence with the spectrum of the fermion theory. However, it is easy to 
see that this correspondence breaks down at high enough energies, if the number 
of fermions N is finite, though it may be very large. In fact, the partition function 
of the fermion theory agrees with that of the collective theory only if we cut-off the 
oscillator frequency of the collective field 15 at N. This cut-off has to be imposed 
by hand; it is not a part of the standard quantization of the collective theory. In 
contrast, as we have argued earlier, in the noncommutative (i.e. exact) formulation 
of the bosonized theory, the spectrum exactly matches with the fermi theory, for 
any number of fermions (including one!) without the need for any cut-off. 

Another interesting quantity to compute is the Wigner density u(p, 8, t) in some 
state. It is possible to write down a manifestly gauge- invariant classical expression 
for u(p,6,t) in terms of the functions p(cr,t) and 0(a,t) which characterize the 
boundary. We have, 

u (p, 6,t) = J da d a 6(a, t) Q(p(a, t) - p) 5{6{a, t)-6). (4.8) 

Here is the familiar step-function. One can easily verify from this expression 
that the equation of motion for u, (d t — de)u = 0, implies the equation of motion 
(|3.8|) for the boundary. To see that this expression satisfies the constraint u 2 = u, 



15 I would like to thank S. Minwalla for pointing this out to me. The calculation of the 
partition function of TV fermions in a harmonic oscillator potential can be equivalently done 
as the calculation of the partition function of a matrix valued harmonic oscillator, gauged under 
U(N). (In fact, the 1/2-BPS sector of TV = 4 super Yang-Mills is actually a U{N) one-matrix 
quantum mechanical system with a harmonic oscillator potential |?6[ [TJ].) The latter calculation 
must take into account only the gauge-invariant states and the partition function over these is 
given by Zn((3) = n^ =1 (l — e"' 3 ™) -1 . Comparing with the well-known partition function of the 
bosonic system described by the Hamiltonian (4.7), namely Z{[3) = n^° =1 (l — e - ' 3 ' 1 ) -1 , we see 
that agreement requires the cut-off. 
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we substitute in it the gauge-fixed solution ( |3.10|) to the equation of motion. Then, 
the a integral in ( |4.8| ) can be explicitly done. We get, 



u(p,9,t) = e(p(t + e)- P ), (4.9) 

which has the desired form. 

Upon quantization, the boundary fluctuates, so in the quantum theory a 
reasonable definition would be in terms of averages in a given state. We define 

u(p,9,t) =< @(p(t + 6) -p) > . (4.10) 

The quantum average can be worked out exactly in any state since we are dealing 
with a free theory. Here we will restrict ourselves to the simplest case of the ground 
state, |0 >. In this case u = u(p) will turn out to be a function of p only. Taking 
a derivative of (|4.10|) with respect to p, using the fourier representation of the 



resulting 5-function, and doing the average gives the result 

d p u{p) = --^ e -^ 2 / 2c , (4.11) 

V 27TC 

where 

oo 

c=< 0\(p(t + 9) - p ) 2 \0 >= h 2 rn (4.12) 



m=l 



To get a finite value for c, one needs to impose a cut-off on the frequency sum by 
hand. Then, using the boundary condition that u(p) vanishes at infinity, we can 
integrate (|4.11|) to get 

u{p) = -^L= r dx e -(*-p°) 2 /2c_ (4. 13 ) 
V2vrc Jp 

As a consistency condition, the density obtained above must satisfy the fixed area 
constraint / dp u(p) = p . This constraint 16 translates into the condition 

Po/VTc = \ r dp (1 - erf(p)), (4.14) 

where erf (p) is the error function. It turns out that this condition is always satisfied 
if iVo < Po/h, where iVo is a cut-off (assumed large) on the sum in ( |4.12|) . Since 
Po/h is finite, though it may be large for small h, we once again see the need to 
impose a cut-off on the oscillator frequencies. 



16 Onc might think that this constraint should be automatically satisfied since, at least formally, 
one has J °° dp u — J °° dp < 0\Q(p(t + 9) — p)\0 >=< 0\p(t + 0)\0 >= po. However, it turns out 
this formal argument does not work since large negative fluctuations can destroy the positivity 
of p{u). Hence the need to impose this constraint explicitly. 
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The classical density function for the ground state is the step-function 6(p — 
p). Quantum corrections work in the right direction and soften the fall-off to an 
exponential. However, for finite N the detailed functional form of the density in 
( |4.13| ) does not match with that of the exact answer for the fermi vacuum given 
in ( |2.16| ). For example, the integrand in the latter case shows rapid oscillations 
with a wavelength which decreases as N~ 1 / 2 for large N . In contrast, the collective 
theory answer for the integrand in ( 4.13Q is a simple gaussian which has no such 
feature. It is important to emphasize here that ( 4.13|) is the exact answer for the 
density in the ground state in the collective theory. There are no corrections. We 
conclude that the collective quantum theory calculation does not agree with the 
exact answer for finite N. 



4.2 Boundary profiles with radially directed tangents 

Fig. 2 (a) shows an example of a boundary profile which has two such points, none 
of which is a point of inflexion. Fig. 2(b) shows an example of a boundary profile 
which has just a point of inflexion. As we have discussed, a gauge cannot be 





(a) 



(b) 



Figure 2: (a) A fluid boundary profile with radially directed tangents at two points, 
(b) A fluid boundary profile with a radially directed tangent at one point which is also 
a point of inflexion. Both droplets are centered at the origin O of the phase plane. 



completely fixed on 9 for such boundaries; gauge-fixing has to be partly done on p. 
As a consequence of this, p does not provide a complete description of the gauge- 
invariant degrees of freedom for such boundaries. In this case, a more convenient 
set of variables is provided by the gauge-invariant quantities a m which are defined 
below: 

1 



a. 



2nh 



da d a 4>(a, t) e 



-im(t+6(cr,t)) 



m 



0,±1,±2, 



(4.15) 



They are manifestly invariant under time-dependent boundary reparametrizations. 
Moreover, using ( |3.5| ) and the gauge-fixed solution to the equations of motion, 
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( CT) , we get 



Oim = ^rl da dj(a) p(a) e~ imS ^ . (4.16) 
27m j 



For 6(a) = a, these are identical to the modes used in the previous subsection 
(with ha = Pq and a_ m = a^) for boundary profiles without any radially directed 
tangents. For the more general boundary profiles under discussion here, there is 
a subtlety in inverting ( [4. 16 ) to express p in terms of these modes. We have from 



(EUD 

TiY^otm e imS{a) = [ da' d a ,d{a') p(a') 5{9{a') - 9(a)). (4.17) 

Let us denote the location of zeroes of d a 9 by ex,, where i = 1, 2, • • ■. For a ^ a i: , 
we get 

fc$> m e mS ^ = sign(«9 >J(a)) p(a), a ^ a t . (4.18) 

m 

However, for a = Oi the right hand side vanishes, leading to the constraint 

]>> m e ime >*)=0, (4.19) 



one for each point ov Equations (|4.18|) and (|4.19|) express one set of gauge-invariant 



variables, namely {{p(a), a ^ (ij}, {^(cTj)}} 17 in terms of another, namely {a m }. 
This latter set clearly provides a more convenient starting point for a gauge- 
invariant description of the quantum dynamics of generic boundary profiles. 

As in the previous subsection, quantization begins with the canonical 
commutation relation ( |4.2| ), which is valid for generic boundary profiles. From 
this we deduce 18 the standard harmonic oscillator commutation relations ( |4.5|) for 
a m . Moreover, as before the Hamiltonian is given by ( [4.61 ). Since p(<Xt) does not 
contribute to the right hand side, we may use ( (4.18| ) in it. This gives precisely 



the expression ( |4.7j ) for the Hamiltonian in terms of the modes. We see that the 
spectrum remains unchanged, and, as before, for finite N it does not agree with 
the spectrum in the fermionic theory. 

Analogous to the computation following ( |4.8|) , here also one may wish to 
compute the Wigner density u(p, 9, t) in some state. First note that (|4.8| ) correctly 
reproduces the classical density for a generic boundary profile. Using the gauge- 
fixed solution ( |3.10| ) to the equation of motion, we get 



u(p, 9,t) = J da d a 9(a) &(p(a) - p) 5(9(a) -t-9). (4.20) 



17 



The complimentary set, namely {{p(ai)}, {0(<j), o ^ 0"^}}, has been gauged away. 
3 One needs to use the relation i / da dj(a) e 4 ( m -") e » = S mn . 
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Let <7fc(t + 9), k = 1,2, ■■• be the points at which the 5-function clicks, i.e. 
9(&k) — t + 9. Then, we get 



u(p,9,t) 



^sign(<9, 



9(a k )) Q(p(a k )-p). 



(4.21) 



The sum excludes points a k = &i at which d a 9 vanishes. It is easy to see that 
the right hand side precisely equals one in the filled region and is zero otherwise. 
The exclusion of zeroes of d a 9 from the sum as well as the presence of the "sign" 
function in this formula is crucial to reproduce the correct answer. 

In the quantum case, we replace the step-function on the right hand side 



of ( 4.21 ) by the average in some state. For the ground state, the average is 
independent of o k and has the value given in ( f4.13|) . The sum is also trivially 
done since the number of solutions, excluding the points where d a 9 vanishes, is 
always odd with one more plus "sign" than the "minus" sign. The net result for 
the density is +1 times that in (|4.13|) , i.e. it is identical to that answer. So as 
there, we conclude that the collective quantum theory calculation does not agree 
with the exact answer for finite N. 

We end this section with the following comment. In the above we have only 
considered configurations with a single droplet centered around the origin in the 
phase plane. For more general configurations consisting of several disconnected 
droplets shown in Fig. 1(b), the density u and hence the action can be written as a 
sum on the different droplets. The origin of phase plane will be inside at most one 
of the droplets. Figs. 3 (a) and (b) show examples of profiles when the origin of the 
phase plane is outside the filled region. In such cases there are always at least two 





(a) 



(b) 



Figure 3: (a) A fluid boundary profile with just the two bounding radially directed 
tangents, (b) A fluid boundary profile with radially directed tangents other than the two 
bounding tangents. The origin O of the phase plane is outside both droplets. 

radially directed tangents. Furthermore, 9 takes a maximum value which is < 2tt. 
By suitably modifying the above discussion to incorporate these differences, one 
can easily extend the present analysis to a generic configuration of droplets. 
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5. Summary and concluding remarks 



The noncommutative theory of the Wigner phase space density developed in 



2£j, |30], [Hp provides an exact bosonization of 2-dimensional non-relativistic 
fermions. Utilizing the construction of noncommutative solitons of [33|, in this 
paper we have shown that the spectrum of the bosonized theory is identical to 
the spectrum of the fermion theory. Moreover, given a Wigner density function 
which solves the equation of motion and constraints of the bosonized theory, we 
have a precise algorithm for building the occupied levels in the corresponding fermi 
state. In contrast, we have shown that the collective quantization of fluctuations 
of fermi fluid droplet boundaries neither reproduces the spectrum nor the details 
of the phase space density, except perhaps for strictly infinite number of fermions. 
Although we can explicitly demonstrate this exactly only for the harmonic oscillator 
potential, we believe the result is valid more generally. 

One possible way out of this disagreement with the collective theory could be 
that the Jacobian arising from the complicated change of variables modifies the 
"classical" action ( |3.6|) at finite N, possibly to all orders in 1/N (or equivalently 
in h). A systematic procedure for arriving at such a "corrected" action would be 
useful, but we have not attempted this here. Even if one is able to find such an 
action, for finite number of fermions one would still need to impose by hand a 
cut-off on the collective field oscillator frequencies at N, which is required for the 
spectrum of the collective theory to match with the spectrum of N fermions. In 
fact, since it is the oscillators defined in ( J4.16|) that are the appropriate gauge- 
invariant variables to use for boundary profiles of any shape, it might be useful to 
reformulate the collective theory directly in terms of N such oscillators, namely 
the set {(a m , a^J, m = 1, 2, • • • , N}. Any "corrected" collective action in terms of 
this set of variables, which reproduces the results of the fermi theory, at least order- 
by-order in H, is likely to have much of the structure of the exact noncommutative 
bosonization. In fact, the latter may be a good staring point for investigating this 
possibility, which we leave for future work. 

The present work has implications for probing quantum gravity in the 1/2- 
BPS sector using supergravity fluctuations and D-branes As discussed 
above, exact equivalence to free fermions via AdS/CFT requires the appearance 
of a noncommutative structure on the gravity side. It would be very interesting if 
one could relate this noncommutative structure to the cut-off required in quantum 
gravity. 
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